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Critical behavior of a traffic flow model 
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The Nagel-Schreckenberg traffic flow model shows a transition from a free flow regime to a jammed 
regime for increasing car density. The measurement of the dynamical structure factor offers the 
chance to observe the evolution of jams without the necessity to define a car to be jammed or 
not. Above the jamming transition the dynamical structure factor exhibits for a given fc-value two 
maxima corresponding to the separation of the system into the free flow phase and jammed phase. 
We obtain from a finite-size scaling analysis of the smallest jam mode that approaching the transition 
long range correlations of the jams occur. 
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I. INTRODUCTION 

Real traffic displays with increasing car density a tran- 
sition from free flow traffic to congested traffic. This 
break down of free flow traffic is accompanied by the oc- 
currence of traffic jams which are not attributable to any 
external cause but only to distance fluctuations between 
following vehicles within very dense and unstable traffic 
(see for instance jij). These traffic jams find expression 
in shock waves, i.e., in backward moving density fluctu- 
ations. This property of jams was already found during 
the 60's in traffic observations 0]. Beyond this early 
experimental investigations recently performed measure- 
ments of real highway traffic leads to the conjecture that 
jams are characterized by several independent parame- 
ters, e.g., the jam velocity, mean flux, etc. gj. 

Since the seminal work of Lighthill and Whitham in 
the middle of the 50's Q many attempts have been made 
to construct more and more sophisticated models which 
incorporate various phenomena occurring in real traffic 
(for an overview see ^Q). Few years ago Nagel and 
Schreckenberg j]]] introduced a cellular automata model, 
which simulates single-lane one-way traffic, and which is 
able to reproduce the main features of traffic flow, back- 
ward moving shock waves and the so-called fundamental 
diagram (see || and references therein). Investigations 
of the Nagel-Schreckenberg traffic flow model show that 
crossing a critical point a transition takes place from a ho- 
mogeneous regime (free flow phase) to an inhomogeneous 
regime which is characterized by a coexistence of the free 
flow phase and the jammed phase PJl0|. Thereby, the 
free flow phase is characterized by a low local density 
and the jammed phase by a high local density, respec- 
tively. Due to the particle conservation of the model 
the transition is realized by the system separating into a 
low density region and a high density region |j] . Despite 
several attempts it is still an open question wether the 
transition can be described as a critical phenomenon (see 
for instance §,[l]-@). 

An attempt to investigate the occurrence of jam was 
made by Vilar and Souza |i~H ] who introduced an order 
parameter which equals the number of standing cars. 
Without noise this quantity vanishes continuously be- 
low the transition. But including noise the amount of 
standing cars is finite for any nonzero density, i.e., the 
average number of standing cars does not vanish below 



the critical density. The crucial point is that below the 
transition a particle may stop due to the noisy slowing 
down rule of the dynamics, but this behavior does not 
coincide with backward moving density fluctuations and 
algebraic decaying correlation functions which would in- 
dicate a critical behavior of the system. At present, a 
convincing definition of an order parameter, which tends 
below the critical density to zero, is not known. Again, it 
is even not clear wether the Nagel-Schreckenberg traffic 
flow model displays criticality at all. 

We use the dynamical structure factor both to make 
predictions about the properties of the different phases 
observed in the Nagel-Schreckenberg model and to inves- 
tigate the question wether the occurrence of jam can be 
described as a phase transition or not. The dynamical 
structure factor which is closely related to the correla- 
tion function is an appropriate tool to do this because 
it naturally distinguishes between the two phases char- 
acterized by positive and negative velocities (free traffic 
flow and backward moving shock waves). Thus the ad- 
vantage of our method of analyses is that the properties 
of both phases can be examined without the necessity to 
define a car to be jammed or not. The paper is organized 
as follows: In Sec. II we briefly describe our analysis of 
the dynamical structure factor and recall the main results 
which were published recently Jl5| . In Sec. Ill we extend 
this method of analysis to the so-called velocity-particle 
space and address the question whether the transition 
from the free flow regime to the phase coexisting regime 
(where the system separates into the jammed and free 
flow phase) can be considered as a phase transition. Bas- 
ing on the dynamical structure factor, our results suggest 
that a continuous phase transition takes place. 



II. MODEL AND SIMULATIONS 

The Nagel-Schreckenberg traffic flow model is based 
on a one-dimensional cellular automaton of linear size L 
and N particles. Integer values describing the position 
r n G {1,2, ...,£}, the velocity v n G {0, 1, w max } and 
the gap g n to the forward neighbor are associated with 
each particle. For each particle, the following update 
steps representing the acceleration, the slowing down, the 
noise, and the motion of the particles are done in parallel: 
(1) if v n < g n the velocity is increased with respect to the 
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maximal velocity, v n — > Mia{v n + l,i> m ax}, (2) to avoid 
crashes the velocity is decreased by v n — > g n if v n > g n , 
(3) if v n > the velocity is decreased by v n — » u„ — 1 
with probability P in order to allow fluctuations, and 
finally (4) the motion of the cars is given by r n — > r n + 
v n . Thus the behavior of the model is determined by 
three parameters, the maximal velocity u ma x, the noise 
parameter P and the global density of cars p = N/L, 
where N denotes the total number of cars and L the 
system size, which is chosen to be L = 32768 throughout 
the whole paper. 

Our analyses is based on the dynamical structure fac- 
tor whose definition is as follows: Consider the occupa- 
tion function 

J 1 if cell r is occupied at time t 
^ ~ \ otherwise. ( > 

The evolution of r] r j leads directly to the space-time di- 
agram where the propagation of theparticles can be vi- 
sualized (see for instance Fig. 2 in [jl6[ ). The dynamical 
structure factor S(k,uj) is then given by 



S(k,cu) = 




where the Fourier transform is taken over a finite rect- 
angle of the space-time diagram of size I x T, i.e., r and 
t are integers ranging from 1 to I and 1 to T, respec- 
tively. Then k and u> are also discrete values k = 2i:m,k jl 
and lo = 2irm UJ /T with 6 {0,1,2,...,/ — 1} and 
m u e {0, 1, 2, T — 1}, respectively. The dynamical 
structure factor S(k, uS) is related to the Fourier trans- 
form of the real space density-density correlation func- 
tion C(r, t) and compared to the analysis of the steady 
state structure factor |J and the related steady state cor- 
relation function j^] it contains both the spatial and the 
temporal evolution of the system. Figure [j] shows the dy- 
namical structure factor both below and above the tran- 
sition. Below the transition S(k,u>) exhibits one mode 
formed by the ridges. This mode is characterized by a 
positive slope v — duj/dk corresponding to the positive 
velocity Vf of the particles in the free flow phase. In- 
creasing the global density, a second mode appears at 
the transition to the coexistence regime. This second 
mode exhibits a negative slope (vj < 0) indicating that it 
corresponds to the backward moving density fluctuations 
in the jammed phase. 

Due to the sign of their characteristic velocities Vf and 
Vj, both phases can clearly be distinguished. Recently 
performed investigations turned out that the character- 
istic velocity of the free flow phase equals the velocity 
(uf = u ma x — P) of free flowing cars in the low density 
limit (p -* 0), i.e., cars of the free flow phase behave as 
independent particles for all densities |lj|. The velocity 
of the jams neither depends on the global density nor on 
the maximum velocity, i.e., it is a function of the noise pa- 
rameter only [^5) . This result coincides with those of re- 
cently published investigations fl7j] , which were obtained 
from a variation analysis of a multi-point-autocorrelation 
function. The jam velocity is, beside the maximum flow, 



the fluxes inside and outside of the jam and the aver- 
age vehicle speeds (see for instance P,|i~8| and references 
therein) a characteristic parameter of real traffic and its 
knowledge is therefore necessary to calibrate any traffic 
model to the conditions observed in real traffic. 



III. THE JAMMED PHASE 

In the last section we briefly described that the anal- 
ysis of the dynamical structure factor allows to deter- 
mine the characteristic velocities of the jammed and the 
free flow phase. In this section we are interested in the 
correlations occurring within the jammed phase. There- 
fore it is necessary to analyse the jammed mode. Unfor- 
tunately, the superposition of the jammed and the free 
flow mode makes it difficult to consider the pure jammed 
mode (see Fig. [j]). Thus a refinement of the analysis 
of the dynamical structure that allows to separate both 
modes is needed. This can be achieved by changing from 
the occupation function rj r< t , defined in real space, to the 
velocity-particle space where the evolution of each parti- 
cle velocity v n ,t is considered. A snapshot of the velocity- 
particle space for a given time is shown in Fig. ^. In 
the free flow phase, where the cars can be considered 
as independent particles ||, the velocities fluctuate ac- 
cording to the noisy slowing down rule between the two 
values w m ax and i) ma x — 1, respectively. Extended regions 
with small or even zero velocities corresponds to traffic 
jams. Comparable to the space-time diagram these re- 
gions move backward in time. The Fourier transform of 
the velocity-particle diagram in two dimensions leads to 



S(k,u) 




FIG. 1. The dynamical structure factor S(k, u) below (up- 
per figure) and above (lower figure) the jamming transition. 
The ridges with maximal S(k,u) indicate the various modes. 
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FIG. 2. Snapshot of the velocity v n ,t as a function of car 
number n at a certain time t and P = 0.5. Jams are char- 
acterized by several particles with low velocity. The free flow 
particles fluctuate independently between v n ,t = W — 1 and 



the dynamical structure factor S v (k,u>) which is given by 

), (3) 



1 



S v (k,Lu) = 



with k = 2nmk/N, to — 2?:m u /T and where the m's 
being integers. 

Compared to the dynamical structure factor of the 
ordinary space-time diagram the dynamical structure 
factor of the velocity-particle space has the advantage 
that the free flow phase contributes only white noise 
(S v (k, ^freeflow = const), i.e., the analysis of the oc- 
curring traffic jams is made easier. Figure || shows the 
structure factor S v (k,u>) above the transition. The struc- 
ture factor displays one ridge with a negative slope, cor- 
responding to the backward moving jams. The notch 
parallel to the /c-axis through oj — is caused by finite- 
size effects and disappears with increasing system size. 
The peak in the center of the diagram [S v (k = 0, oj = 0)] 
describes the velocity fluctuations of the whole system. 
It is an average over both coexisting phases and therefore 
it leads to uncertain results on the occurrence of jams. 

In the following we carry out a quantitative analysis 
of the ridges to obtain information about correlations 
between jammed particles. We will show that among 
those particles long range correlations exist only above 



S„(*,w) 




FIG. 3. The dynamical structure factor S v (k, oj) of the ve- 
locity-particle space above the transition for P — 0.5. Notice 
that the free flow phase contributes to S v (k,u}) only white 



the transition whereas below the transition these corre- 
lations are restricted to a finite correlation length. In 
Fig. |^ we plot the values of the dynamical structure fac- 
tor S v (k,oj) for uj/k = Vj, i.e., the values along the ridge 
which correspond to the jam modes. In order to take 
finite-size effects into account, we use the scaling ansatz 



S v (k,u>)\ 



ut/k- 



N- e f(N c k). 



(4) 



For 9 = —1.22 and £ = 0.34 we obtain a convincing data 
collapse of all curves. The dynamical structure factor 
decays algebraically, 



(5) 



where the exponent is given by 7 w 3.6 ± 0.2. This al- 
gebraic decay of the dynamical structure factor indicates 
that the corresponding correlation function is also char- 
acterized by an algebraic decay, i.e., the system displays 
long range correlations above the transition. Next we are 
interested in the values of S v (k, oj) for u)/k = Vj below the 
critical density. These values are shown in Fig. 0. Due to 
the finite amount of standing cars below the transition, 
the jam mode does not vanish. Our analysis shows that 
the jam mode decreases like a Lorentz curve, 



S v{k,0j)\ w/k=v . 



1 



i + (fc0 2 



(6) 



Herein £ denotes a correlation length, defined in the 
velocity-particle space. The term c takes into consider- 
ation that, caused by the free flow phase, the jam mode 
does not tend to zero for large k but to a finite value. 
Beside the jam mode, Fig. || shows a Lorentz curve ac- 
cording to Eq. (||), which has been fitted to the data. 
To estimate the influence of finite-size effects on the cor- 
relation length, we first determined the £ for different 
values of N (with T — N) and found that £ is almost 
in-affected by finite-size effects for the densities shown in 
Fig. 0. Therefore, these values of the correlation length 
equal within negligible errors those values of the corre- 
lation length found in the limit N — > 00 (with N = T). 
The dependence of the correlation length on p c — p is 
shown in the inset of Fig. |f| With p approaching p c the 
correlation length diverges as 



£ ~ (p c -p)-", 



(7) 



with v = 0.92 ±0.05. 

Since £ describes the correlation of particles of the 
jammed phase, the occurrence of long range order in the 
system at p = p c indicates that the system displays crit- 
icality. As mentioned above, the correlation length £ is 
infinite above the transition (long range order) and finite 
below. This coincides with measurements of the life time 
distribution of jams at the critical point fl9|| which dis- 
plays a power-law behavior, i.e., traffic jams occur on all 
time scales. Below the critical density the life time of the 
jams is finite and no long range correlations of jams can 
occur. 

In the following we address the question how long range 
order occurs by approaching the transition point. There- 
fore we consider the smallest positive mode on the ridge 
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FIG. 4. The scaling plot of S v (k, u>) along the jam mode for 
various values of TV (with TV £ {128, 256, 512, 1024, 2048} and 
TV = T) and a fixed value of the global density with p > p c . 
For 6 — —1.22 and £ = 0.34 we got a good data collapse 
indicating that the correlation of the jams decay algebraically 
(dashed line with slope 7 = —3.6). 
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FIG. 5. The dynamical structure factor S v (k,u>) for 
uj/k = Vj, -Umax = 5, P = 0.5 and TV = 1024. The dashed 
line corresponds to a Lorentzian fit according to Eq. . The 
inset shows the correlation length £ as a function of p c —p. The 
obtained power-law behavior (solid line) agrees with Eq. (Q). 



of jammed particles, S v (k — * 0,cj — > 0) with w/fc = Uj. 
This value contains the information how long range and 
long time correlations appear when the transition takes 
place and it is believed to be closely related to an or- 
der parameter (see for instance pof). If the transition 
from the free flow regime to the jammed regime can be 
described as a critical phenomenon, S v (k — > 0, u> — > 0) 
should vanish below the transition (p < p c ). Simulating 
finite system sizes this means that the smallest jam mode 
obeys the finite-size scaling ansatz 

S v (k -» 0, lj -» 0) = (TVT)-^ 2 /[(7VT)*/ 2 (p - Pc )], (8) 

with w/fc = Uj. We measured SV^fc — > 0,cj — > 0) for 
various values of TV and T and for P = 0.321 which cor- 
responds to a jam velocity vj rj —1/2 The finite-size 
scaling works and for y w —0.2 and a; « 0.1 we got a good 
data collapse, which is plotted in Fig. [| Especially we 
obtain the data collapse for different ratios of TV and T, 
which shows that the system behaves isotropic in TV and 
T. To investigate the dependence of the exponents on P, 
we determined them also for another value of the noise 
parameter (P = 0.519 corresponding to Vj « —1/3 ). 
Due to the size of the corresponding error bars no signif- 
icant P dependence of the exponents could be observed. 
Further investigation with improved accuracy are needed 
to clarify this point. 

Note that it is not justified to identify the scaling expo- 
nents [Eq. (||)] with the usual critical exponents of second 
order phase transitions (y = fi/v and x = 1/f). First it is 
not clear if S v (k — > 0, ui — > 0) equals the order parameter 
(see and references therein) . The second point is that 
the usual finite-size scaling ansatz rests on the validity of 
the hyperscaling relation between the exponents (see for 
instance [pl| ). Despite these restrictions the finite-size 
scaling analysis of the smallest jam mode reveals that it 



vanishes below the transition. In the hydrodynamic limit 
(N — > co, T — > 00) no long range correlations in space 
and time occur for p < p c . Above the critical value the 
correlation function displays an algebraic decay. Since 
the correlations of the jams are finite below and infinite 
above p c , the system displays critical behavior. This 
agrees with the above mentioned investigations of the life 
time distribution [flfjf of jams and with measurements of 
the relaxation time of the system. It was shown that 
the relaxation time r diverges at the critical point with 
increasing system size r ~ L z with a P depending ex- 
ponent z [[12|-|l4|. But one has to mention that above 
the transition the measurements of the relaxation time 
yields unphysical result, in the sense that the relaxation 
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(P-Pc) ^ 

FIG. 6. The smallest mode of S v {k -> 0,lu -> 0), with 
u/k = Vj, for different values of TV (TV G {128, 256, 512, 1024}) 
and ratios TV/T (TV/T G {0.5, 1, 2, 4}) and P = 0.321. 
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time becomes negative |14| . We think that the origin of 
this behavior is caused by the inhomogeneous character 
of the system above the transition where the system sep- 
arates into two coexisting phases. The relaxation time 
measurement does not take this inhomogeneous charac- 
ter into account. 



IV. CONCLUSIONS 

We studied numerically the Nagel-Schreckenberg traf- 
fic flow model. The investigation of the dynamical struc- 
ture factor allowed us to examine the transition of the 
system from a free flow regime to a jammed regime. 
Above the transition the dynamical structure factor ex- 
hibits two modes corresponding to the coexisting free flow 
and jammed phase. Due to the sign of their characteris- 
tic velocities Vf and Vj , both phases can clearly be distin- 
guished. 

The analysis of the dynamical structure factor of the 
velocity-particle space shows that above the transition 
the system exhibits long range correlations of the jammed 
particles and below the transition the correlations display 
a finite correlation length that diverges at the critical 
point, indicating that a continuous phase transition takes 
place. Using a finite-size scaling analysis we showed that 
in the hydrodynamic limit the smallest jam mode which 
corresponds to the long range correlations of jams van- 
ishes below the transition. We think that an extended 
investigation of this quantity could lead to a convinc- 
ing definition of an order parameter which describes the 
transition from the free flow to the jammed regime. 
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